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A regular effective approach in studying the many- 
particle problem at the microscopic level is the quantum- 
field approach borrowed by the statistical physics of non- 
relativistic systems from quantum field theory. This 
method was employed by Matsubara fH to construct 
the thermodynamic perturbation theory and the diagram 
technique for investigating many-particle systems at fi- 
nite temperatures. The Matsubara's approach was im- 
proved by Abrikosov, Gor'kov, and Dzyaloshinskii 0, Q 
and Fradkin Q . The technique developed in these works 
is applicable for systems, whose symmetry coincides with 
the symmetry of their Hamiltonians. In order to de- 
scribe the states of many-particle systems with sponta- 
neously broken symmetry (superfluidity, superconductiv- 
ity, ferromagnetism, crystals, etc.) and the phase tran- 
sitions between the states with various symmetries, the 
quantum-field approach is to be modified. In spite of 
the progress achieved in the development of the theory 
of many-particle systems that are in the state with spon- 
taneously broken symmetry, the further development of 
the ways to describe such systems on the basis of fun- 
damental physical principles is a challenging direction in 
the evolution of quantum statistical physics. 

The approach proposed in this work is based on a con- 
sistent account, already in the zero-order approximation, 
of the essential features of the system state under con- 
sideration; first of all, its symmetry. It can be made con- 
sistently only in the case where the interaction between 
particles is made allowance for, at least to some extent, 
not later than in the main approximation. A requirement 
for the initial approximation to be simple dictates the ne- 
cessity to choose the SCF model as the main approxima- 
tion. The relative simplicity of this model is defined by 
the circumstance that it preserves, as well as the model 
of independent particles, the one-particle (to say more 
accurate, quasi-one-particle) description of the system as 
a gas composed of quasiparticles which are characterized 
by individual wave functions. An important element of 
the theory is making use of the Bogolyubov concept of 
quasiaverages 0, [i'] which has been formulated, taking 



into account the self-consistent approach, in work Q. 

The SCF model is known to be widely used for calcu- 
lating atomic spectra jS], the structure of atomic nuclei 
9] , and the properties of molecules and solids ; there- 
fore, one may hope that the properties of many-particle 
systems can be calculated well enough by this method, 
with the quality of such calculations growing as computer 
facilities are being developed. The main attention in this 
work is given to the development, on the basis of the 
generalized SCF model, of the quantum-field approach 
to describe many-particle Fermi systems which are in 
states with spontaneously broken symmetry (including 
spatially nonhomogeneous ones) and to the construction 
of a perturbation theory and a diagram technique. It 
has been shown that the thermodynamic functions found 
from a microscopic consideration in the framework of the 
SCF model obey the correct thermodynamic relation- 
ships, and, therefore, the SCF model, which is funda- 
mental for the theory proposed, is not self-contradictory. 
It should be emphasized that the approach, which is pro- 
posed for the consideration of many-particle systems, is 
based only upon the general principles of non-relativistic 
quantum theory and statistical mechanics and does not 
require additional hypotheses. 

1 . The motion of a single particle in the external field 
Uq{y) is described by the Schrodinger equation 



dq'Ho{q,q')ipj{q') = efVjC?) , 



(1) 



where the notation q = {r, a} stands for the set of spa- 
tial, r, and discrete spin, cr, coordinates; the integration 
includes the summation over spin indices; the subscript 
j comprises the full set of quantum numbers, including 
the spin projection, which characterizes the state of an 
individual particle; (pj (q) is the wave function of the par- 
ticle; and e^°^ is its energy. The kernel in Eq. ([T]) looks 
like 



Ho{q,q')(pj{q') 



2to 



AS{q-q') + Uo{r)S{q-q'), 



where S{q — q') = S{r — r') S^a', m is the particle's mass. 
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and A is the Laplacian. We will use the methods of 
secondary quantization, having defined the operators of 
creation, a^, and annihilation, Uj, of real Fermi particles. 
Let us also define the field operators 

= "f^iq) =J2^*AlH ■ (2) 



The particle operators and field operators ([2]) obey the 
notorious anticommutation relations 6j. The Hamilto- 
nian of a many-particle Fermi system, expressed in terms 
of the field operators, looks like 



H = Jdqdq'^+{q)Hiq,q')^{q') + 

+ 1 I dqdq' ^+{q)^+{q')U{r,r')^{q')^{q) . 



(3) 



where 



H{q,q')^Ho{q,q')-fiS{q-q'). 



and [/(r, r') is the potential of two-particle interaction 
independent of spin variables. While studying many- 
particle systems with broken symmetry, it is convenient 
to suppose the considered system to be in contact with 
a thermostat and to have an opportunity to exchange 
both energy and particles with it; i.e. the total energy 
and the total number of particles in the system are not 
regarded fixed. The thermostat is characterized by two 
parameters: the temperature T and the chemical poten- 
tial n- In the state of thermodynamic equilibrium, the 
same parameters characterize the system of particles as 
well. Thus, we use the grand canonical ensemble and 
deal with the Hamiltonian that includes the term with 
the chemical potential — /iiV, where N is the operator of 
the number of particles. 

2. To pass to the SCF model, let us present the initial 
Hamiltonian ([3|) as a sum of two terms: 



H — Ho + He , 



(4) 



where the first term is the Hamiltonian of the SCF model 
which is quadratic in the field operators, 

H„ = jdq dq'l^^+{q)[H{q, q') + W{q, q')^{q')] + 

+ ivl/+(5)A(g,g')*+('z') + ^*(9')A*(g, + ^^0, 

(5) 

and the second one is the correlation Hamiltonian 
He = dqdq'l^^+{q)^+{q')U{r,r')^{q')^iq)- 
-2^+{q)Wiq,q')^{q')- 



-^+{q)A{q,q')M> + iq') - ^ (q' ) A* {q , q' )^ {q)\ ~ E',, 



which takes those interactions into consideration that 
were not accounted for in the SCF approximation. For- 
mulae (HI and dH) include self-consistent, still unknown 
potentials W{q, q') and A(q, q') which satisfy, owing to 
the self-consistency of the Hamiltonian, the conditions 

W{q,q')^W*iq\q), A{q, q') = -Aiq' , q), (7) 

as well as the non-operator term Eq , the choice of which is 
essential for the thermodynamics of the considered model 
to be constructed correctly. Therefore, in the SCF model, 
Hamiltonian ([3]) is replaced by the simpler model Hamil- 
tonian Hq The latter contains potentials which will 
be determined from the condition of the best approxima- 
tion of the exact Hamiltonian H by the model Hamil- 
tonian Hq. An important qualitative difference between 
these two Hamiltonians consists in that the initial Hamil- 
tonian does not depend on the system state, whereas 
the self-consistent Hamiltonian Hq, as will be shown, de- 
pends on the system state and the thermodynamic vari- 
ables through the self-consistent potentials W{q,q') and 
A{q,q'). It is this property of the self-consistent Hamil- 
tonian that allows one to describe the states with broken 
symmetry. When constructing the perturbation theory 
for many-particle systems with broken symmetry, it is 
natural that the self-consistent Hamiltonian Hq should 
be chosen as the basic one, and the correlation Hamilto- 
nian He as a perturbation. 

With the help of Bogolyubov's canonical u-v- 
transformations 



i 
i 

we can diagonalize the self-consistent Hamiltonian (0) 



(8) 



Hq — Eq 



It li ■ 



(9) 



(6) 



Here, Eq is the non-operator part of the Hamiltonian; 
Ei is the energy of elementary excitations, quasiparticles, 
reckoned from the chemical potential level; i is the full 
set of quantum numbers, including the spin projection, 
which characterizes the quasiparticle state; and the oper- 
ators "ff and 7i describe the processes of creation and an- 
nihilation of quasiparticles. The quasiparticle description 
is widely used in condensed mater physics. In the SCF 
model, the notion of quasiparticles, which possess the in- 
finite lifetime in this approximation, arises naturally as a 
result of the reduction of Hamiltonian ([5]) to the diago- 
nal form ([9]). The relative simplicity of such a model lies 
in that it retains the one-particle (to say more precisely, 
one-quasiparticle) description of the system. The set of 
factors {ui{q), Vi{q)} is the two-component wave function 
of a quasiparticle. For the transition from the initial self- 
consistent Hamiltonian ([5]) to the diagonalized one (O to 
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be feasible, the factors in canonical transformations ^ 
must satisfy the Bogolyubov-de Gennes system of equa- 
tions [H, 12| which looks, in the most general case, like 



dq' [n{q, q')uM') + A(q, q')v,{q')] = U,{q) , (10) 



dq' [n*iq,q')v,{q') + A* {q, q')u,{q')] = -e,v,{q) , 

(11) 

where Vl{q,q') = H{q,q') + W{q,q'). The requirement 
that for transformations (|S]) to be canonical results in the 
normalization and completeness of the solutions of the 
self-consistent equations (fTOj) and ([TT|) [13]. The average 
values of operators in the SCF model are expressed in 
terms of the normal, p, and abnormal, r, one-particle 
density matrices: 

p{q,q') = {'f+{q')^{q))o = 
= ^[u,(<z)<(q') /, + v:{q)v,{q'){l /,)] , (12) 

i 

T{q,q') = {^{q')^{q))o = 



E 



u^{q)v:{q')h+vUq)u.{q'){l-fd], ^^^^ 



where the quasiparticle distribution function fi has the 
same form as in the model of ideal gas, 



= = (exp/3£, + l)-i 



(14) 



and /3 — 1/T is the inverse temperature. Since the quasi- 
particle energy is the functional of fi, formula (|14p turns 
out a complicated nonlinear equation for the distribution 
function, being similar to that which takes place in the 
Landau phcnomcnological theory of a Fermi liquid [3]. 
In formulae (fT2|) and (fT3|) . the averaging is carried on 
with the statistical operator 



(15) 



po = cxp^(rio - Hq) , 

where the normalizing constant 51o = — TlnSpe 
determined by the condition Sp po = 1 and has the mean- 
ing of the thermodynamic potential of the system in the 
SCF model. 

In order that the system of equations (fTOl) and ([TT|l be 
completely defined, one has to express the self-consistent 
potentials ([7]) in terms of the functions u{q) and v{q). 
Using the variational principle [13], we find the connec- 
tion between the self-consistent potentials with the one- 
particle density matrices (|12p and ([T3| : 



(16) 



Wiq,q') = ~U{r,r')piq,q') + 
+ 5{q-q') [dq"U{v.v")p{q'\q") 



Substituting Eqs. (HH) and ^ into Eqs. ^ and ^ 
gives rise to a closed system of nonlinear integro- 
differential equations for the wave functions u{q) and 
v{q): 



U{r) - p. ~ ei + 



+ dq'U{v,v')piq\q') 



Ut{q)' 



dq' U{v, v')[p{q, q')u,{q') - r{q, q')v^{q')] = , 



(18) 



U{y) -p + ei + 



+ dq'U{v,v')p{q\q') 



■"i{q)- 



dq' C/(r, r')[p*(g, q')v,{q') - T*{q, q')uM')] = . 



(19) 

The chemical potential p is related to the average number 
of particles N as 



Jd^rn{r), n(r) = ^ p((z, g) 



(20) 



where n(r) is the particle number density. The system of 
equations (IT8| and p9)) , together with relations and 
(pn| , describes fermion systems at finite temperatures in 
the SCF approximation, being also valid when the sym- 
metry of the system states is spontaneously broken. In 
particular, these equations are applicable for describing 
both magnetic properties and superfluid (superconduct- 
ing for charged particles) states. For spatially nonhomo- 
geneous states, where the characteristic variation length 
of the functions u(q) and v{q) is much longer than the ef- 
fective range of the interparticle potential, Eqs. and 
(IT9l) can be presented in the differential form. The con- 
stant E'q in Eq. ^ is defined by the formula [l^] 

E'o^-l fdqdq'U{r,r'){^+{q)^+{q')^{q')^{q))o = 



\ jdq dq'U{r,r')x 



A{q,q')^U{r,r')r{q,q'). 



(17) 



|P(9,<?')P-P(9,9W,9')-W9,<?')P 

(21) 

In the SCF approximation, the average value of the exact 
Hamiltonian is equal to the average value of the self- 
consistent Hamiltonian: {H)q = {Hq)o. 

In many cases, there is no necessity in knowing the 
quasiparticle wave functions in order to find the equilib- 
rium characteristics of the researched system; the one- 
particle density matrices will do. From Eqs. (1181) and 
(IT^ . as well as formulae and p^ . the system of 
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equations for the one-particle density matrices reads: 



ing Eq. (|24l) into account, we obtain 



— [Ap(g,q') - AV(g,<z')] - Pir) - U{v'Mq,qr 



dq" U{vy)-U{v'y) p{q,q")p{q",q') 



-pW, q")p{q, q) + ^(g, q")T* {q", q) 



0, 



(22) 



2m 



[AT(g,q') + A'T(g,g')]- 
[U{v) + U{v') + t/(r, r') - 2y] riq, q') + 

dq" \u{r, r") + U(r\ r")] \p{q, q")T{q" , q'y 



-Piq",q")riq,q')+Tiq,q")piq',q")\ =0. 

(23) 

The solutions of this system and, hence, the correspond- 
ing self-consistent fields (fT6|) and ([T7|) can possess differ- 
ent symmetries, including that which is lower than the 
symmetry of the initial Hamiltonian ^ , and thus can de- 
scribe states with a spontaneously broken symmetry. In 
particular, the system of equations ([^^ and (1^51) has both 
normal solutions, for which T{q,q') = and p{q,q') 0, 
and "superfluid" ones, for which both r(g, q') and p{q, q') 
do not vanish. The anomalous density matrix T{q,q') or 
the self-consistent potential A{q, q') can be considered as 
the microscopic order parameters of the superfluid state. 
If the dependencies of the one-particle density matrices 
on the spin variables are such that p{q,q') oc Saa' and 
T{q,q') oc cr^^l (a^^^ is the Pauli spin matrix), the sys- 
tem is invariant with respect to spin rotations. Other- 
wise, there is a magnetic ordering in the many-particle 
system. If the thermodynamic parameters are fixed, only 
that state among the possible ones of the system will be 
realized really, whose thermodynamic potential is mini- 
mal. 

3. A distinctive feature of the SCF model, which 
should be taken into account when deriving thermo- 
dynamic relations from Hamiltonian ([S]), is that this 
Hamiltonian contains self-consistent potentials and an 
operator-free term which depend on the temperature 
and the chemical potential. Only the correct choice of 
the self-consistent potentials W and A and the quantity 
Eq ensures that the thermodynamic relations would be 
satisfied. Using the definitions of thermodynamic po- 
tential and entropy 5*0 = — Sp(polnpo), one can 
demonstrate [isl . [iBj that the thermodynamic relation 
Hq — E — TSq — fj,N, where E is the total energy of 
the system, is satisfied, and the variation of the thermo- 
dynamic potential is equal to the averaged variation of 



Srto = {SHo)o 



(24) 



By varying the self-consistent Hamiltonian which is ex- 
pressed in terms of one-particle density matrices and tak- 



SHo 



5^0 



Sp{q,q') \5p{q,q')/ Q ST*{q,q') \(5r* (g, q')/ o ^' 

(25) 

The relations (|16p and p7)) between the potentials W 
and A and the one-particle density matrices, which have 
been established with the help of the variational princi- 
ple, make the thermodynamic potential extremal, as is 
seen from Eq. (j25p . with respect to its variation over the 
one-particle density matrices p and r. Due to Eq. (|25p . 
in the case where the volume of the system is fixed, the 
ordinary thermodynamic relation 



dflo 



-SodT - Ndp 



(26) 



is satisfied. The total energy can be found either by aver- 
aging the energy operator directly or with the help of the 
thermodynamic relation in terms of the thermodynamic 
potential 



(27) 



According to Eqs. (|26l) and ^7}, the fact that the self- 
consistent Hamiltonian involves the potentials depending 
on thermodynamic variables does not violate thermody- 
namic relations, as might have appeared [l6| . so that 
the SCF approximation in statistics is intrinsically non- 
contradictory. 

The total energy of the system of Fermi particles in the 
SCF model is a sum of several contributions (the kinetic 
energy K, energy of particles in an external field Ue, 
energy of direct particle- particle interaction Ud, energy 
of exchange interaction Ucx, and energy of condensation 
into a superfluid state Uc) and can be written down in 
the form 

i i 

(28) 

One can readily see that the total energy is not the sum 
of energies of individual quasiparticles. Carrying out av- 
eraging in Eq. lU and taking into account Eq. ([SS)) , we 
obtain the constant Eq in the diagonalized Hamiltonian 
®: 

E = -{UD + U,y,+Uc) - Y^e, jdq\v,{q)?. (29) 

i 

Now, the ultimate form of the thermodynamic potential 
in the SCF approximation can be found easily as 



f^O = -(C/d + C/cx + C/c)-^£. Jdq\v,{q)\^- 



(30) 
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In formulae IMl - ISO 



Ud^^J d^r Sr' C/(r, r') 7i(r) n(r') , 



dqdq'U{r,r')\p{q,q')\\ 



Uc^- Jdqdq'U{r,r')\Tiq,q')\\ 

The total number of quasiparticles is defined by the 
relation 

^. = E/- (31) 

i 

while the total number of particles looks like 

N = N, + J2th^ fdq\vM\'- (32) 
i •' 

Hence, in the normal (non-superfluid) state the number 
of particles is equal to the number of quasiparticles, sim- 
ilarly to what takes place in the theory of a normal Lan- 
dau Fermi liquid In the superfluid state, as follows 
from Eq. (P^ . the number of quasiparticles is always less 
than the number of particles. It can be interpreted in a 
way that a certain number of particles in the superfluid 
phase forms the condensate of Cooper pairs and does not 
contribute to the formation of quasiparticle excitations. 

The relations of the SCF theory quoted above can be 
presented 17| in the form of the relations of the Fermi 
liquid theory [l^. The equations of the generalized SCF 
model formulated here lead to the results of the BCS 
theory of superconductivity [3], as well as the rela- 
tions of the theory of a superfluid Fermi liquid which 
has been constructed in works 1^ 2^ . The equations of 
the generalized SCF model are capable to describe mag- 
netic properties of collective electrons and bring about 

Note that the 
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the relations of the Stoner theory 
proposed theory makes it possible to write down the 
known Stoner criterion for ferromagnetism in another 
form, Ip{so) > 1, where / is the parameter of exchange 
interaction and /c(eo) is the density of states at the Fermi 
level. In the case of point-like interaction between parti- 
cles, C/(r,r') = UoS{r — r'), this criterion looks like 



I 



(33) 



where I = {V/NY^^ is the average distance between par- 
ticles, a — Uom/{4TTti^) is the scattering length, and 

— ^{^y^^ ~ 0.5 is a numerical parameter. The ap- 
proach proposed here enables the competition of super- 
conductivity and magnetism, as well as spatially modu- 
lated states, to be studied. 

4. While studying theoretically the systems with bro- 
ken symmetry, one can not take advantage of the conven- 
tional definition of the average for calculating observable 



characteristics, because the symmetry of a state with bro- 
ken symmetry is lower than that of its Hamiltonian. At 
the same time, when calculating the averages according 
to the routines of statistical mechanics, the symmetry 
of the averages coincides with that of the Hamiltonian. 
Such a contradiction does not arise in the SCF model, 
because the system of self-consistent equations has solu- 
tions with symmetry lower than that of the initial Hamil- 
tonian. To overcome those difficulties, Bogolyubov Q 
introduced the concept of quasiaverages into statistical 
mechanics. According to this concept, the averages in 
states with broken symmetry should be calculated not 
using Hamiltonian ([3]) but a Hamiltonian, which differs 
from ([3]) by terms that violate its symmetry in an appro- 
priate way. In the framework of such an approach, how- 
ever, some uncertainty in the introduction of fields that 
violate symmetry remains. As the choice of symmetry- 
violating fields does not depend on interparticle interac- 
tions, it may happen that the interactions do not allow 
the existence of states possessing the symmetry imposed 
by the field introduced. A way to determine quasiaver- 
ages using the self-consistent Hamiltonian as an additive 
that violates the symmetry was proposed in work [7,]. In 
this case, the system can possess only the symmetry that 
is allowed by interparticle interactions. 

Though the symmetry of each of the Hamiltonians Hq 
and He, which depend on the system state, can be lower 
than the symmetry of the initial Hamiltonian, the sym- 
metry of H does not depend, naturally, on the way how 
it was split and, thus, remains invariant. Therefore, in 
order to describe the systems with broken symmetry, let 
us introduce a more general Hamiltonian 



Hg — Ho + gHc 



(34) 



which depends on a real positive parameter g. It is ob- 
vious that, at g = 1, this Hamiltonian coincides with 
the initial one ^ and, at 5 = 0, with the self-consistent 
Hamiltonian ([5]). The variation of this parameter from 
zero to unity means the inclusion of the correlation in- 
teraction. If g is very close to unity, Hamiltonian (|34|) 
practically coincides with the initial one ([3]). The ma- 
jor difference, however, consists in that the symmetry of 
Hamiltonian (p4)) coincides at that with the symmetry 
of the self-consistent Hamiltonian and can be lower than 
the symmetry of the initial Hamiltonian. Let us define 
the statistical operator 



Pg 



-Tin Spe 



(35) 



The quasiaverage of an arbitrary operator will be deter- 
mined by the relation 



(A) = lim lim Sp pgA . 



(36) 



If the values of the thermodynamic variables p and T are 
fixed, quasiaverages (j36l) can be not equal to conventional 
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averages and, therefore, cannot describe the states with 
broken symmetry. From the mathematical point of view, 
a possible divergence between averages and quasiaverages 
is known to arise from the dependence of the result 

on the sequence of the passages to the limit in Eq. ([5^ . 
The passage to the limit of the "coupling constant" g 
must be carried out after the thermodynamic passage to 
the limits V ^ oo and N ^ oo, provided N/V = const. 
If the symmetry is not broken, quasiaverages p6p are 
identical to the relevant conventional averages. 

5. Correlation Hamiltonian ^ which we have selected 
as a perturbation has a rather complicated structure. 
However, it can be written down in a much more compact 
form, and the perturbation theory will accept a simpler 
form if one uses the concept of the normal product of op- 
erators. This concept plays an essential role in quantum 
field theory |,23i] . In the temperature-involved technique 
which was put forward in works [11-0] j the concept of 
the normal product is not applied; therefore, the analogy 
with quantum field theory is incomplete. 

To define the normal product of Fermi operators, one 
has to pass preliminarily to the particle-hole representa- 
tion. Such a transition at T = does not meet difficulties 
[l^ . At non-zero temperatures, this procedure is not so 
obvious, because one-particle states are not divided un- 
ambiguously into occupied and free ones; so that every 
state may be either free or occupied with a certain prob- 
ability. However, the concept of the normal product can 
be generalized for finite temperatures so that it will be 
independent of the transition to the particle-hole repre- 
sentation. 

For the further consideration, it is convenient to in- 
troduce the notation for operators, using the "isotopic" 
index which acquires two values, 1 and 2: 

(37) 

We also use the notation a which means 

J 1, if a ^2, , - 

"=|2, ifa^l. (3^) 

The temperature normal product of two operators is 
defined by the formula 

iV(*i*2) = *1*2 -P21, (39) 

where p2i = (^'1^2)0 is the one-particle density matrix 
and 1 = (qi, ii). It is evident that 

(iV(^'i^'2))o =0. (40) 

Let us give the general definition of the normal product 
of operators which should be valid for both the Fermi 
and Bose statistics. We introduce the notion of operator 
pairing which means the average over a self-consistent 
state, 

viV2 = {ViV2)o , (41) 



where iji is any of the operators aai,'^ai or jai- The 
product of an arbitrary number of operators containing 
the pairings is defined as 

= 0(771 ?73)o(?7fcr7m)oX (42) 
X ?72?74 . . . rjk^ir-jk+i ■ ■ ■ rjjn-iVm+i ■ ■ ■ Vj-iVj j 

where a is a numerical factor which is equal to unity for 
Bose operators and (—1)^ for Fermi ones; and p is the 
number of permutations that are needed for putting the 
paired operators side by side in the initial order. Taking 
into account the given definition of pairings, the normal 
product of an arbitrary number of operators is defined 
by the formula 

N{r]iri2 . .-Vj) = mV2 ■ --Vj- 

-v1v2m ■ ■ - m - vtmvl ■ ■ - m- 

— (all other products with one pairing)-!- 

+ (all other products with two pairings) — . . . . 

(43) 

Thus, the temperature normal product of operators is 
defined as a sum of operator products, which include all 
possible pairings (including the term without pairings). 
The sigh plus is selected if the number of pairings in the 
product is even, and the minus if odd. The average of the 
A'^-product of any number of operators in the Schrodinger 
or interaction representation over a self-consistent state 
is equal to zero, 

(7V(*i...*,)o = 0, (44) 

except the average of the A^-product of a c-number which 
is N{c) = c by definition. We do not develop the general 
proof of property (j44|) in detail, but it is easy to check 
in a straightforward manner that it is fulfilled for the 
A^-product of, e.g., four operators: 

-*1*2P43 - *3*4/521 - *2^'3/541 - *l*4/532 + 
+ \I'2*4/031 + *1*3P42 - P31P42 + P21PAZ + P41P32- 

(45) 

An important property of the SCF model is that it al- 
lows a rather complicated correlation Hamiltonian ^ to 
be represented as the normal product of four field oper- 
ators: 

Hc = \ jdqdq' U{vy) N[^+{q)^+{q')^{q')^{q)\ = 

= i E U{v, v')N[^^{q)^J^q')^d.q')'^d.q)\ ■ 

(46) 



7 



Similarly, this Hamiltonian can be expressed in terms of 
quasiparticle operators: 



t^i234 A^(7i 727374) , 



(47) 



1234 



where the matrix element of the interaction potential 



C/i234 = Jdqdq' f/(r,r') u^^^qK:' iQ^CiqK:' (q) 

(48) 

is expressed through the coefficients of the Bogolyubov 
transformation matrix 



Vt{q) u*{q) 



Thus, the Hamiltonian of the correlation interaction can 
be written down in the form of a normal product not 
only for the normal systems at zero temperature [l6| but 
also for the states with spontaneously broken symmetry 
at finite temperatures. 

6. An L-point temperature Green's function (GF) is 
defined as 



G{l,2,...L)^i'^{TrA{l)A{2)...A{L) 



(49) 



where the averaging should be understood as the quasi- 
averaging p6p . and each number stands for the whole set 
of variables. The operators that are averaged in Eq. 
are in the Heisenberg-Matsubara representation: 



(50) 



where Aa is the corresponding operator in the 
Schrodinger representation, r the Matsubara "time" pa- 
rameter (0 < T < /3), and Tt the operator of chronolog- 
ical ordering P, d. If i(l) = ^-(1), Eq. (03) defines an 
i-point field GF, and if A{1) =7(1), an L-point quasi- 
particle GF. When studying Fermi systems, only GFs 
with even L's are considered. In particular, a two-point 
(one-particle) GFs are defined by the formulae 



G""' (gr, q'r') = -{T, i> ^{qr) i> ^, {q'r')) , 
G^^'{it,i't') = -(T,7m(T)7a'»'(T')) ■ 



(51) 



These functions are 4x4 matrices in the spin and "iso- 
topic" spaces. The components of GFs ([5T|) that are 
diagonal in isotopic indices are anomalous and distinct 
from zero only in the superfluid state, while the non- 
diagonal ones are non-zero both in the superfluid and 
normal states. To construct the perturbation theory, the 
operators in the Matsubara interaction representation are 
to be introduced: 



-rHa 



(52) 



Using these operators, we define the temperature GFs in 
the framework of the SCF model by relations 



G(")""V,9V) = -(r,M'„(qr)vI/„,(gV))„, 
G(")""V,*V) = -(T,7«(r)7o'.'(T'))o- 



(53) 



Here, the averaging is carried on with the statistical op- 
erator (|15p over the self-consistent state. Functions ([5T|) 
and (15^ depend only on the "time" difference r — r' and 
satisfy the symmetry relations 



G""(<z,g';r-r') = -G""(g',<z;r'-r) = 
= G"" *(g', g; T - r') = -G""'*(g, g'; r' - r) 



(54) 



For GFs ([55]) in the SCF model, the closed equations 
are valid. It is most convenient to calculate the quasi- 
particle GF, for which 



5G(0)""'(ir, i'r') 



dr 



r;„e,G(")""'(zT,^V) 



(55) 



where rja = (—!)"■ Expanding the GF in Eq. ([55l) in a 
Fourier series, we find the Fourier component 



^ii' y^nj — ; 

lUJn + Va£i 



(56) 



where u;„ =7rT(2n 4- 1). Therefore, both in the normal 
and superfluid states, only those terms are distinct from 
zero which are non-diagonal in isotopic indices. Taking 
into account Eq. (jS]), we find the GF in the field repre- 
sentation: 

i.a" 

In this case, distinct from zero are the terms which are 
either diagonal or non-diagonal in isotopic indices. 

7. Let us construct the diagram technique to find the 
field and quasiparticle GFs. Define the contraction of 
operators in the interaction representation by the relation 



n {q, r)n' (q, r') = T.(vl/„(g, r)*^, (g', t'))- 
-N{^^(q,T)-^^,{q',T')). 



(58) 



Contraction (|58|) is a c-number which coincides, with an 
accuracy to a sign, with the GF 

■^liq, t)-^1, {q\ t') = -G(°)""'(qr, q'r'). (59) 



We note that product p3|) and contraction ((58)) are so 
defined that Wick's theorems remain valid in that form as 
they were formulated in quantum field theory ; there- 
fore, the perturbation theory can be built in the standard 
way 0, [3 H^l ■ To present the diagram technique in a 
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more compact form, we designate variables, on which the 
GF depends, by a single number, e.g., 1 = (ri, ai, ci, ri) 
and I = (ri, ai, (Ti, Ti), and consider the integration over 
the digital variable as the integration over all continuous 
coordinates and the summation over all discrete ones. In 
order to construct the perturbation theory, it is necessary 
to express the temperature GFs in terms of operators in 
the interaction representation and to pass from the av- 
eraging with the statistical operator pg [see Eq. ([35]) ] to 
that with the statistical operator po [see Eq. ([T5|) ]. Car- 
rying out such a transition [l6| , we obtain 



(60) 



for the one-particle GF, where the temperature scattering 
matrix looks like 



o'(/3) — Tr exp 



r0 

/ Hc{T')dT' 
JO 



(61) 



where Hc{T)~e'^^° He e^"^^" is the perturbation Hamil- 
tonian in the interaction representation. Since the per- 
turbation Hamiltonians and (|T7|) are presented in 
the normally ordered form, the averaging over the self- 
consistent state is to be applied to the T-products of op- 
erators that are combined in fours under the symbol of 
the A^-product. This reduces the number of pairings sub- 
stantially and, accordingly, simplifies the diagram tech- 
nique. In principle, one could construct the theory with- 
out resorting to the normal ordering of operators. In this 
case, the theory would contain a large number of "su- 
perfluous" diagrams which would not contribute to the 
final result, being reduced in every order of perturbation 
theory. Using the A^-ordered form of the perturbation 
Hamiltonian allows the appearance of such diagrams to 
be excluded. In this case, similarly to the standard tech- 
nique [l,[lB: 24 1, the theorem on connectivity is valid, so 
that the denominator in Eq. (pO]) should not be taken into 
account, and only connected diagrams should be allowed 
for in the numerator. 

We note that the total thermodynamic potential of the 
system is expressed in terms of the temperature scatter- 
ing matrix averaged over the self-consistent state. This 
averaged quantity can be presented in the form [ij, [l^l 



{<7{l3))o = exp 



E^'^»(/3))c 



.n=0 



SO that the total thermodynamic potential is defined by 
the formula 

oc 

ri = O0-rE((T„(/3))0con„. (62) 

Expanding the exponent in Eq. (|6ip in a series, an arbi- 
trary GF can be written down as 



G(l,2,...L) = £G(")(l,2,...i) 



(63) 



The contributions of the n-th order to the thermody- 
namic potential and the GF are determined by the rela- 
tions 



o"f-l)" 

(cT„(/?))oconn^ ^ ^ , ' / . 



(64) 



dTn{Tr Hc{ti) . . . i?c(Tn))oc 



G<")(1,2,...£)^ ■'^"^'V"'' ^ 



dT[... 



dT,;(T, A(1)A(2) . ..A{L)Hc{t[) . . .i?C«))oconn. 

(65) 

Owing to the properties of the normal product, 
(f i(/3))o = 0; therefore, the contribution of corrections to 
the thermodynamic potential in the SCF approximation 
arises only in the second order of perturbation, i.e. the 
summation in formula ()62|) starts from n = 2. Anal- 
ogously, the first-order contribution to the total one- 
particle GF equals zero, so that the correction to the GF 
calculated in the framework of the SCF model emerges 
only in the second order of perturbation theory. 

First, let us formulate the diagram technique for the 
field GFs. We introduce the following graphic notations: 

1 « 2 for G(°\i,2) = -^'''(1)*"(2), 

1 2 for ;7(l,2) = ;7(ri,r2)5(Ti-r2)x 

^(^aiQ2 ^0102)' 

Since a fermionic GF is antisymmetric with respect to the 
index permutation, G^^\l,2) — —G^^\2,l), in order to 
take this circumstance into account, we define the direc- 
tion of the Green's line and agree that the index pointed 
by the arrow stands first in the analytical record of the 
GF. The line, whose arrow is directed in the opposite 
side, corresponds to to the same GF but with the minus 
sign. Each vertex of the interaction line (dashed line) 
joins two Green's lines. Since Hamiltonian includes, 
in pairs, operators that differ in only their isotopic in- 
dices (e.g., vl/(l') and *(1')), those two GFs include the 
vertex indices which differ in a bar above. 

Let us formulate the rules for calculating the contribu- 
tion of the n-th order to the L-point temperature GF: 

1) draw 2n vertices linked in pairs by dashed lines; 

2) link all the vertices with solid Green's lines in all topo- 
logically nonequivalent ways, so that one Green's line 
should enter into each vertex of the interaction line and 
one should leave it; 

3) L/2 vertices must be the sinks for the external Green's 
lines, while the other L/2 vertices must be the sources for 
them; every incoming external Green's line is linked by 
a sequence of solid lines to one of outgoing Green's lines; 
each such solid line is directed away from the external 



9 



a) 



2" 



\2' 



b) 2 



^2 



Figure 1. 



1" 



1^ 




Figure 2. 



<— 2 



index, the position of which in the analytical record of 
the GF is on the right; 

4) the vertices that are not the sinks or the sources for 
the external Green's lines are linked by closed solid lines 
with arbitrary directions; 

5) no Green's line can link vertices which belong to the 
same dashed line; 

6) confront the graphic representations with their analyt- 
ical expressions: each solid line with the GF G'^°\l,2), 
and each dashed line with the interaction potential 
;^(1,2); 

7) the first index of the GF, which was confronted with 
a solid line that enters into the vertex in question, must 
be taken with a bar; 

8) carry out the integration and summation over all the 
variables that were put in accordance with each vertex; 

9) the analytical expression that has been constructed in 
the indicated way according to the diagram should be 
multiplied by the factor 



n+k+C 



where n is the order of the diagram; k the total number 
of loops; ^2 the number of loops which pass through two 
vertices of interaction lines; and Q the number of permu- 
tations needed for the external indices linked by a solid 
line to be arranged in the order as they enter into the 
analytical record of the GF. 

It is important to emphasize that, in this technique, no 
Green's line can link vertices belonging to the same inter- 
action line (item 5). Such elements are known 
to be involved into the diagram technique based on the 
ideal gas approximation. Those diagrams can include 
both Green's lines leaving and entering the same ver- 
tex ("bubbles") and Green's lines lin king the vertices of 
the same interaction line ( "oysters" ) 25| . The diagrams 



which contain such elements are impossible in the tech- 
nique proposed, because, owing to the normal form of 
the perturbation Hamiltonian, there are no contractions 
between operators that stand under the sign of the same 
iV-product. The absence of such diagrams in this tech- 
nique is also natural, because the diagrams made up of 
these elements define the SCF approximation which has 
already been taken into account in this case in the main 
approximation. In Fig.[l] as an example, the diagrams of 
the second order are shown which define the corrections 
to the one-particle GF in the field representation. 



The formulated diagram technique, as well the tech- 
nique based on the approximation of non-interacting par- 
ticles '^l admits the diagrams to be summed up by sepa- 
rate blocks and the graphic methods of summation to be 
used. 

The diagrams for the temperature scattering matrix 
are constructed by the same rules, as those for the con- 
struction of GFs. The former differ from the latter by the 
absence of external lines. As is known 0, [l^], it does not 
allow the graphic summation of the infinite sequences of 
diagrams to be carried out in this case. 

For the practical use of the diagram technique, it is 
more convenient to pass to the frequency representation. 
In this case, the rules of the diagram technique undergo 
the following modifications: 

1) every Green's line is associated with Fourier- 
component ([57l) G^^\l, 2;a;„), and every external incom- 
ing line should be associated with a frequency with the 
minus sign; 

2) every dashed line is associated with the potential 

^(1,2) = ;7(ri,r2)((5aiQ2 - 

3) the frequency conservation law must be fulfilled: the 
sum of the frequencies of Green's lines which enter the 
end points of every dashed line of the interaction is equal 
to the sum of frequencies of the outgoing lines, which is 
taken into account by introducing the multiplier A(w„j^ -|- 
a;„2 — — tiJm)', with A(a;) = 1 if cj = 0, and A(w) = 
otherwise; 

4) the additional multiplier emerges before the 
expression that correspond to the diagram. 

Now, let us formulate the rules of the diagram tech- 
nique in the case where the quasiparticle description is 
used. Taking into account that the adjacent operators 
can be permutated under the sign of normal product with 
changing the sign, Hamiltonian (|T7l) can be represented 
in the form 



He 



4! 



^^(71727374) 



(66) 



1234 



where the antisymmetrized potential 

t^l234 = ^^1234 + ^^1342 + C^1423 + t^2314 + £^2431 + ^^3412 " 
— ^^1243 ~ ^^1324 — U 1432 — U2d.ii — C/2413 — Uzi21- 



(67) 
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was introduced. Let us introduce a notation for the po- 
tential 



Uil, 2, 3, 4) = Ui234 - ^2) <5(ti - T3) 5[t^ - t^) , 

(68) 

which is antisymmetric with respect to the permutation 
of variables. A single digit in the potential C/(l,2,3,4) 
designates the following set of variables: 1 = (zi, ai, ti). 
Every square in the diagram is confronted with the ma- 
trix element of the interaction potential 



C/(l,2,3,4) 



while every Green's line with the zero-order approxima- 
tion of the GF in the quasiparticle representation. Indices 
at the square vertices must be arranged either clockwise 
or counterclockwise in that order as they enter into the 
expression for the matrix element (j68p . 

Thus, the rules for calculating the contribution of the 
n-th order to the Green's quasiparticle function are as 
follows: 

1) draw n squares that correspond to the matrix elements 
of the interaction potential; 

2) link the vertices of different squares with solid lines 
in all such topologically nonequivalents ways, that two 
vertices of a square serve as sinks for one Green's line 
each and two other vertices of the same square be the 
sources of one Green's line each; 

3) confront the graphic representations with their ana- 
lytical expressions; sum up over the variables related to 
each vertex and integrate over "time" variables; 

4) put the multiplier (— l)"+'^+'3/P„ before the expres- 
sion that has been constructed according to the diagram, 
where n is the order of the diagram, k the number of 
closed loops in it, Q the number of permutations which 
are necessary to arrange the external indices linked by a 
solid line in that order as they enter into the analytical 
record of the GF, P„ the number of permutations of the 
indices in the square vertices that do not result in new 
expressions. 

In the quasiparticle representation, we get only one 
diagram of the second order which is shown in Fig. [2] 

For the frequency representation of quasiparticle GFs, 
the rules are modified in the same way as in the case of 
field GFs. 

8. Similarly as in the quantum-field approach that uses 
the model of independent particles as the zero-order ap- 
proximation, the concept of self-energy and vertex parts 
can be introduced in the approach that is developed here, 
and the Dyson's equation that couples those two func- 
tions can be derived. The equation for a one-particle GF 



can be presented in the following form: 
5G(1,2) 



5ti 



-5(1-2)- 



- j d3 [H{1, 3) + W{1, 3) + S(l, 3)] G(3, 2) , 



(69) 



where the self-energy function I](l,2) is defined by the 
relation 

d3 S(l, 3) G(3, 2)^^Jd3 U{1, 3) x 

X [G(l, 3, 3, 2) + G^"\3, 3) G(l, 2) -2 G(°\3, 1) G(3, 2)] . 

(70) 

Formula (l69l) includes the functions 



(71) 



i?(l,3) ^ Ha^aAHi 13)^(^1 ~ T3) , 

M^(l,3) = M/a,a3(gi,g3)5(ri-T3), 



where 



Waa'{q,q') = 



H{q,q') 
-H*{q,q') 



A{q,q') W{q,q') 
-W*{q,q') -A*{q,q') 



With regard for the equations for a GF in the SCF ap- 
proximation, we obtain the known relation for the self- 
energy function: 

G(l, 2) = G("\i, 2) + / d3d4 G(l, 3) E(3, 4) G(°\4, 2) . 



(72) 

The vertex part r(l, 2, 3, 4) is defined by the formula 
G(l,2,3,4) = 
= G(l, 2) G(3, 4) -I- G(l, 4) G(2, 3) - G(l, 3) G(2, A)+ 

+ y"dl'd2'd3'd4'r(l',2',3',4')x 

X G(r, 1) G(2', 2) G(3', 3) G(4', 4) . 

(73) 

From Eq. (1701) . taking Eq. ((75)) into account, we ob- 
tain the Dyson's equation which establishes a relation 
between the self-energy and vertex parts in Fermi sys- 
tems with spontaneously broken symmetry: 

E(l, 2) = i 1^3 U{1, 3){[G(3, 3) - GW(3, 3)] S{1 - 2)- 

-2[G(3,l)-G(°\3,l)]5(2-3)} + 
1 



(74) 



^ J dl'd2'd3'd3U{l,3)x 
xr(l', 2', 3', 2) G(i', 1) G(2', 3) G(3', 3) , 
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The poles of the vertex function introduced by relation 
((73)1 define the dispersion law of collective excitations in 
the many-body system. 

The methods of quantum field theory applied in sta- 
tistical physics were extended in this work to describe 
any eligible states in non-relativistic Fermi systems with 
spontaneously broken symmetry. We managed to do this, 
mainly, owing to two circumstances. First, it is the SCF 
model formulated in the most general form that was used 
as the main approximation. Secondly, it is the procedure 
for calculating the quasiaverages which uses the fields 
that are defined by this model. It is essential that, in 
the given approach, the correlation Hamiltonian consid- 
ered as a perturbation can be presented in the normal 
form, which allows a plenty of diagrams not contribut- 
ing to the final result to be excluded from consideration 
and the diagram technique to be presented in a compact 
form. The approach suggested does not contain any as- 
sumptions and is based only on the general principles 
of quantum mechanics and statistical physics. It can 
be applied for regular researches of equilibrium proper- 
ties of many-particle systems with spontaneously broken 
symmetry (magnetically and spatially ordered, supercon- 
ducting, superfluid, etc. systems) and the phenomena in 
them at a microscopic level. The method proposed can 
be extended onto the description of non-relativistic Bose 
systems with spontaneously broken symmetry, in partic- 
ularjjof superfluid systems with broken phase symmetry 
H, [23. In author's opinion, the general approach de- 
veloped in this work can also be effectively used for a 
proper description of states with spontaneously broken 
symmetry in the relativistic field theory and the theory 
of elementary particles 
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